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Weak asymptotic solution of the phase field system 
in the case of confiuence of free boundaries in the 
Stefan problem with undercooling 

V.G. Danilov* 


Abstract 

We assume that the Stefan problem with undercooling has a classical so¬ 
lution until the moment of contact of free boundaries and the free boundaries 
have continuous velocities until the moment of contact. Under these assump¬ 
tions, we construct a smooth approximation of the global solution of the Stefan 
problem with undercooling, which, until the contact, gives the classical solu¬ 
tion mentioned above and, after the contact, gives a solution which is the 
solution of the heat equation. 


1 Introduction 

In this paper, we study the confluence of free boundaries in the Stefan problem 
with undercooling in the one-dimensional case. We at once note that an 
analysis of the one-dimensional problem is not the final goal but is only a 
necessary step in the study of the multidimensional problem. 

In the preset paper, we study the problem in the domain Q = Q x [0, U], 
where U = [li, I 2 ] is an interval. We assume that the interval U is divided into 
three parts and Q~{t) as follows: 

n~{t) = [tpi{t),(p2{t)], U+(t) = [(/?2(t),^2], 
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where ipi{t), i = 1,2, are the free boundaries of phases ”+” and We 

assume that the phase ”+” occupies the intervals and the phase ” 

occupies the interval 

We shall construct a smooth approximation of solutions of the Stefan 
problem with undercooling under the assumption that the motion of the free 
boundary is the motion of the front of a nonlinear wave and the confluence 
of free boundaries is interaction of solitary nonlinear waves. The possibility 
of this interpretation is given by the models of phase held [Ij proposed by 
G. Caginalp. In fact, the choice of the method for approximating the limit 
Stefan problems with undercooling is unessential for us, because we do not 
prove that the approximations thus constructed are close to the corresponding 
solutions of the phase held system. 

For example, we could use the dehnition of the generalized solution of the 
limit problems including the order function (the nonlinear wave) |l()j . 

Here our considerations are based on the following simple fact. Suppose 
that there are two families of solutions (exact and approximate solutions) of 
some problem depending on a small parameter e. Suppose that both these 
families have the properties that permit passing to the limit as e —> 0 in the 
weak sense. Suppose also that the family of approximate solutions satishes 
a problem with a right-hand side small as e —> 0 in the weak sense. Then 
the weak limits of both families are solutions of the same limit problem and 
if the latter has a unique solution, then the difference between the exact and 
approximate solutions tends to 0 as e ^ 0 at least in the weak sense. 

We recall that the phase held system has the form 


Le = 


du 


u — V? 

eLu - 

e 


-0 = 0 , 


( 1 ) 


where 

L — — — ^^2 ) X G [li, I 2 ], t G [0, ti]. 

The function 6 = 0(x,t,e) has the meaning of the temperature, and the 
function u = u{x,t,e), which is called the order function, determined the 
phase state of the medium: u ~ — 1 corresponds to the phase ” in the 
domain H“(t), and tt ~ 1 corresponds to the phase ”+” in the domains 

Passing to the limit as e —> 0 in o, we obtain the Stefan problem with 
undercooling. 

This passage to the limit is possible, for example, in the case where the 
corresponding limit problems have classical solutions. In this case, the weak 
limits as e —> 0 of solutions m give these solutions laiiiiTsi. 
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The smooth approximations of solutions of the Stefan problem with ki¬ 
netic undercooling constructed in this paper are approximate (in the above 
sense) solutions of the phase field system and they admit a weak passage to 
the limit as e ^ 0. In this case, we obtain the limit problems (and their so¬ 
lutions) describing the process of confluence of free boundaries. To construct 
these approximations, we use the assumption that the classical sharp fronted 
solution of the Stefan problem with kinetic undercooling exists until the con¬ 
fluence of free boundaries begins. This is a natural assumption in our paper 
(otherwise, it is not clear the confluence of what is considered) and can be 
proved^. 

We repeat that, under this assumption, we construct an approximation 
of solution of the limit problem, which is smooth for e > 0 and uniformly 
bounded in e ^ 0. As is well known, in a similar problem about the propa¬ 
gation of shock waves, the existence of such an approximation distinguishes a 
unique solution. 

By t* G [0, fi) we denote the instant of confluence of free boundaries. 
Then, for any t ^ t* — 6 for all d > 0, we see that the limit problems have 
solutions. The asymptotic solution of system o has the form 


qas 


U 


as 

£ 




1 -hcuo 

+ 0JQ 


— X + 




e 7 

X - (P2{t)' 




( 2 ) 

( 3 ) 


Here uJi{z) —> 0,1 as z ^ Too, G §(1^^) for /c > 0, <,3i(f), i = 

1,2, are smooth functions, (pi < (p 2 , ^o{z) = tanh( 2 :), and Lo{t,zi,Z 2 ) G 
C“([0, f*]; §(M^)). By 8(1^*^) we denote the Schwartz space of smooth rapidly 
decreasing functions. If the initial data for o has the form ©, (H at f = 0, 
then, for t ^ t* — 6, we have the estimate 


||u - ur; C’CO, T; l2(]r1)|| t \\e - er; £\Q)\\ ^ ce'^, /i ^ 3/2, 

where {9,u) is a solution of system (1) (see ^E])- Here Q = H x [0,f* — <5), 
and the constant c is independent of e. 

The main obstacle to the construction of solutions of the form m, 
which could be used to describe the confluence of free boundaries, is the fact 

^A. M. Meirmanov. A private communication. 
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that, instead of an ordinary differential equation whose solution is the function 
ujo{z) mis], in the case of confluence of free boundaries, we must deal with a 
partial differential equation for which the explicit form of the exact solution 
is unknown. 

In the present paper, we use the technique of the weak asymptotics method 
mi, which allows us to avoid this problem. Let us explain several basic 
points. 

Definition 1. A family of functions f{x, t, e) G L^{Q) integrable with respect 
to X for all t G [0,ti] and for e > 0 admits the estimate Op/(e^) if, for any 
test function ((x) G C^, we have the estimate 



f{x,t,e)C{x)dx 




( 4 ) 


where the constant depends on ti and the test function C(33)- 

Generalizing 0 , we shall say that the family of distributions f{x,t,e) 
depending on t and e as on parameters admits the estimate C>x)/(e^) if, for 
any test function C{x)i we have the estimate 




Example 1. Let uj{z) G S(]R^), uj{z) dz = 1, and let xq G [li,l2]- Then we 
have 

uj ^ — 6{x) = O'D'ie), oj ^= 6{x — Xq) J uj{z) dz + Ox)'(e)- 


Example 2. Suppose that uJi{z) G C°°, {oji)' G §(M^), lim^^-oo = 0, 
and lim^^+oo = 1, i = 1,2. Then we have 

(a) 

tOi - H{x- Xi) = O-p'is). 

(b) 



= Bi H{x - xi) + B 2 H{x - X 2 ) + C>D'(e), 


where Ax = xi - X 2 , Bi{p) G (7°°, B'-{p) G S(M^), Bi{p) + B 2 {p) = 1, 
Bi{+oo) = 1, and i?i(— 00 ) = 0. 


4 








Example 3. The preceding relations easily imply the formula 


H{x — x\)H{x — X2) = B - J H{x — xi) + (1 — B)H{x — X2) + C>2?'(e), 

where B{p) G C°^ is the function Bi from the preceding example. 

Example 4. The following relation is a corollary of Definition 1: 

( 7 \ m 

Ov'{e^) = Ov'{e^) 

for all m G Z+ and a > 0. Generally speaking, this is not true for the 
derivatives w.r.t. t. 


The relations in Examples 2(a), 2(b) and 4 are obvious, but, for the 
reader’s convenience, we prove some of these formulas in Section 5. 

Here we only note that, in view of items (a) and (b), one can say that 
the product of approximations of the Heaviside functions or of the Heaviside 
functions themselves is a linear combination with accuracy up to small terms 
(~ Ox>'(^))' general nonlinear functions of linear combinations of approx¬ 
imations of the Heaviside functions can be linearized similarly. Precisely this 
property underlies the constructive study of the interaction between nonlinear 
waves with localized fast variations. 


Definition 2. A pair of smooth functions {it, 9) is a weak asymptotic solution 
of the phase field system o if, for any test functions ({x),^{x) G C^{Q), the 
following relations hold; 


[ {ut + 9t)Cdx+ [ O^C^dx = 0{e), 

: / UtUxi dx + ^ I dx , 

Jn ^ Jn ^ J 


u « 1 1 ^ J 

X-Y + jlf'* 


+ 77 / dx - 

./o OX 


^G (0,1/2). 


(5) 

( 6 ) 


The left-hand side of Eq. is obtained by multiplying the second equation 
in system dU) by Ux and integrating by parts. The reminders 0(e“), a = 1,/i, 
in the right-hand sides of © and must be bounded locally in t, i.e., for 
t G [0, ti], we have 


max 10(e)I ^ Ct,e, 


Ct-^ = const. 


This construction was introduced and analyzed in [^. 
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The fact that a noninteger exponent appears in the right-hand side of ® 
is not directly related to the technique of the weak asymptotics method, see 
Examples 1-4. The source of this noninteger exponent is the nonsmoothness 
of the function 0 (the temperature), which appears at the instant of confluence 
of the free boundaries. 

Relations (5)-(6) can be rewritten as 


a - t \ 

^t + ut - + Ojy (e) 




u 


eut = e-^ + u{l - -it) + Ov'{£^)- 


The paper is organized as follows. 

In Section 2, we explain the structure of ansatzes of approximations of the 
temperature and the order function. These ansatzes are constructed under the 
assumption the classical sharp fronted solution of the Stefan problem exists. 
At the beginning of Section 2, we formulate what we exactly need. 

Next, in Section 3, we substitute the constructed ansatzes in system (1) 
and derive equations for the unknown functions contained in the ansatzes. 

The results of these calculations are summarized in Theorems 1-3 in Sec¬ 
tion 3. The final result of this paper is formulated in Theorem 4 in Section 3. 
It states that the assumption on the existence of the classical solution is suffi¬ 
cient for constructing formulas for the weak asymptotic solution of system (1) 
in the sense of Definition 2. Next, in Section 4, we analyze the constructed 
formulas and derive the following effects: 

(a) the weak asymptotic solution is smooth for t > t*, the absolute values 
of the free boundaries velocities are equal to each other at the contact 
moment; 

(b) the temperature has a negative jump at the instant and at the point of 

confluence of the free boundaries, and this jump is equal to the half-sum 
of the limits of the velocities of the free boundaries as t —> — 0. 

In particular, it follows from (a) and (b) that the velocities of the free 
boundaries have jumps at the point of contact. 

These effects can also be discovered in the numerical analysis of the process 
of confluence of free boundaries. Some helpful technical results are given in 
Section 5. We note that the only example known to the authors, where the 
confluence of the free boundaries is studied, is given in HU. 

In general, this paper turned out to be technically rather complicated and 
long in spite of the fact that the details in several justifications in Section 5 
were omitted. 
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Although, as was shown above, the results obtained here are a necessary 
step in the study of the multidimensional problem, this paper shows that the 
following classification can be introduced: 

(1) problems in which the confluence of free boundaries leads to disap¬ 
pearance of one of the phases; 

(2) problems in which the domain occupied by one of the phases changes 
its connectivity, but the number of phases remains the same. 

In this paper, we consider an example precisely from the first class of 
problems. As was mentioned above, we do not justify the asymptotics of the 
constructed solution. But this can be done based on our constructions. It is 
easy to see that the main role here is played by the (not proved) existence of 
the classical solution up to the moment of confluence of the free boundaries. 
This assumption reduces justifying the constructed weak asymptotic solution 
to estimating the soluton of the heat equation with the right-hand side /g 
admitting the estimate 


fe = Oj,:{en, /iG (0,1/2) 

and with zero initial and boundary conditions. An analysis of the structnre 
of this right-hand side shows that fs ■ as e —> 0 is a linear combination 
of functions 5'{x — (fi) and 6{x — ipi), i = 1,2, with coefficients depending on 
t, T, and these coefficients converge fast to zero as r —> ±oo. 

Hence we can conclude that for x ^ the solution of this heat equation 
belongs to C“ for e > 0 and admits the estimate O(e^). In the whole domain 
H X [0,T], a rough analysis based on general theorems |14j shows that the 
solution belongs to d > 0, and has the estimate 0{s^) in the norm of 

this space. In particular, the weak limit of the constructed weak asymptotic 
solution is equal to the exact generalized global solution of the heat equation 
in the phase field system. 


2 Ansatz of the approximation of the so¬ 
lution of the Stefan problem with kinetic 
undercooling 


We recall that, in our case, the Stefan problem with kinetic undercooling has 


the form 


86 _ 8^6 

8t 8x‘^ ’ 


x e [li,l 2 ], ipi, i = 1,2, 


(7) 


7 


(8) 

( 9 ) 


= 2(-l)*+V*t. 

Relations ©-© are supplemented with (Dirichlet or Neumann) boundary 
conditions for x = U, i = 1,2, and with a consistent initial condition. Obvi¬ 
ously, problem 0-® can be written as 

LO = 2if!2t5{x - (p2)-2<f>it6[x - >^i), (10) 

«L.«, = 

We shall assume that the initial conditions are chosen so that the problem in 
question has the classical solution, (^i(O) < <,32(0), and there exists a t = t* G 
[0,ti] such that — 0) = <^ 2 (i* — 0). More precisely, we assume that 

(i) the limits -o exist, i = 1,2; 

(ii) 


'de' 

dx 


9{x,t) e (^[o,r),c2 untit) 

n {c^ {n-{t)) u (sd+(t)) u (o+(t))) n c{n) 


(here D denotes the interior of the domain D), cf. |3]. 

As will be shown below, these assumptions permit describing the interac¬ 
tion (the confluence of free boundaries) constructively. 

We shall construct this approximation (ansatz) as a weak asymptotic so¬ 
lution of system (1). First, we introduce the ansatz of the order function. It 
has the form 



1 + Wo 



+ Wo 


/3 


X-ip2 


- Wo 




( 11 ) 


Here the unknowns are the functions j3 and ipi = <^i(t, e), i = 1, 2. In what 
follows, we write them more precisely. Now we only note that if /? > 0, then, 
for <y9i < ^p 2 i the function ii coincides up to 0{e^) with the sum of the first 
three terms in the right-hand side of ©, for ^pi = <^ 2 , we have it = 1-|-C>x>'(e), 
and for ipi > ip 2 , we have ti = 1 -|- 0{e^). These relations can easily be 
verified directly. We note that, under the above assumptions, we can continue 
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the functions (pi{t) to the interval [0, ti] preserving the smoothness and the 
sign of the derivatives. We choose such continuations and denote them by 
^Pio{t). Now, we can write the functions ipi{t,£) and (3{t,e) more precisely. 
Namely, we set 

(Pi{t,e) = ( 12 ) 

where 

= <P2oit) - y^io{t), T = i)Ql£. ( 13 ) 

We note that 


T ^ OO, t < t*, 

T ^ — OO, t > t*, 

as £ ^ 0. Thus, the value of the variable r characterizes the process of 
confluence of free boundaries. Indeed, we could set r = (t* — t)/e, but, for 
the next formulas become simpler (this can be seen from the formulas 
in Example 2 if we put Xi = (piQ{t)). 

Similarly, we set 

P{t) = 1 +j3i{T). (14) 

In addition, we assume that 

^ r^oo, a = 0,1. (15) 

0 < <5i < Pi{t) < 62 , Si, 62 = const, 

ifii ^ ip~i, ipii-^p~-^ = 0 {\t\~^), t^-oo, = const. 

Assumptions m for the functions (fn are the assumption that as e —> 0 the 
functions (pi{t,e) approximate continuous functions with a possible disconti¬ 
nuity (jump) of the derivatives at t = t*. 

More precisely, the functions ipi{t,£) determined in (12) are approxima¬ 
tions of the functions 


(Pi{t,0) = ipio + ^poiPi^Hi-ipo), 

where = lim-r^-ooipi{t,£) — (/?j(t,0) = 0(e), and H[z) is the 
Heaviside function. 

Now we describe the ansatz of the approximation of the temperature 
9{x,t). Here we also must use the continuation procedure for matching the 
temperature for t < t*, when there are three subdomains of the domain Q, 
and the temperature for t > t*, when there is only one phase. For this, we 
first introduce a “model” of the temperature whose continuation is reduced to 
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the problem of continuation of functions depending only on the time t. This 
model must have the same structure as the temperature 6. Obviously, the 
simplest function of this form is a function linear in x in the domains 
and quadratic in the domain 
Namely, we set 


d9 


dx 


a;=(^2±0 


de 


dx 


x=<pi±0 


±7^ 


We note that, by assumptions (i) and (ii), the functions j = 1,2, 
depend on t, are continuous for t < t*, and have limits as t ^ t*—0. Therefore, 
they can be continued to the interval [0,ti], ti > t*, so that the properties 
of being continuous and the signs are preserved. We shall use the previous 
notation for the continued functions. Now we can introduce the temperature 
model 


^’ = 7 i - x)H{ifi - x) + {x - ip2)H{x - ip2) 

+ ^~{x,t)— —— —H{x-(pi)H{(p2-x) (16) 

+ ^{x,t)— - - x)H{x - (P2) + I, 

where 7“, 7, and I are linear functions of x, 7“!^^^^. = 7~(t), = 7i(t), 

and I\x=tpi ~ III more detail, the functions j{x,t) and 7“(x,t) 

can be written in the form 

- 11-^2 t *^ 7 ^^ + 72 " 

7 

71+72 , *7 71-72 * +1++2 

(fit - ^2t , *\+lt + +2i 

1 =-X — X )-. 

2 2 

We note that, by Lemma 3 in Section 5 , for t ^ t*, we have the relations 

ipi{t,e) = ipi{t) + 0(e), i = l,2. (17) 
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This implies that, for t < t*, 


'df 


'de' 



dx 

X=ipi 

dx 

1 

\x=(pi 

'dT 

dx 

II 

1 

dO' 

dx 

6(x - pi) + OD'{e 


Using the formula in Example 3 above, we rewrite the product 

H{x- ipi)H{ip 2 - x) 


(18) 


as follows 


H{x- (fi)H{(f 2 -x) = B{p)[H{x -(pi) - H{x- P 2 )] + Ov'{e). (19) 

H{(pi - x)H{x - (P 2 ) = (1 - B{p))[H{x - (P 2 ) - H{x - (pi)] + Ojxie). 

Here B{p) G B{p) ^ 1, p ^ 00 , B{p) 0, p ^ - 00 , = 0{\p\-^), 

IpI ^ 00 for any A'" > 0, a > 0. 

In what follows, we write B{t) instead of B{p). The difference is that B{p) 
is an unknonw function (p is unknown), while B{t) is a known function. The 
replacement of B{p) by B{t) in (19) and (20) implies a correction of order 
Ox)'(e) in the right-hand sides. 

But we do not use this remark and, by definition, we set 


^ = 71 i‘fi - x)H{pi - x) 72 (x - ip 2 )H{x - P 2 ) (20) 

+ J~{x,t)— —— —B{t)[H{x - pi) - H{x - P 2 )] 

+ j{x,t)— -—^(1 - B{t))[H{x - P 2 ) - H{x - pi)] +1. 

Now we can write the expression uniformly in tp = p 2 — pi- We have 


5 x 2 

+ & B(rm. - ,0 - (X - ^.)1 

+ ^ (1 - B(rnm. - - H(. - ^,)l 

+ [7^ +12 B{t) + 72(1 - B{t))] 5 {x - P2) 

+ [7r + 'ytB{ t) - 71(1 - ^(t))] 6{x - pi) + Ox,/(e), 
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where Fi (x, t, r) is a bounded piecewise continuous function. We set 7 ^ = —72 
and 7 ^ = 7 i and see that the coefficients for 5{x — (fi) in ( 1211 ) take the form 

[12 + I 2 ] B{t)6{x - if2) + [ 71 “ + it] B{t)5{x - v?i). (22) 

Now we note that the inequalities 

[H{x - ifi) - H{x - ip 2 )\ > 0, [H{x -(P 2 ) - H{x- v?i)] /ip < 0, 

hold for 7 ^ 0. We recall that, by (d and the Stefan condition 

‘fiot = 7i^ + 7r > <^20t = -(7^ + 72") 

for t < t* and assume that the continuations of the functions contained in 
these relations are chosen so that these relations hold for t £ 

The assumption that 6{x, t) is the classical solution of problem 0-® 
implies that 

9-f = qeC{[0,t*y,C\n)) (23) 

and the limit q exists as t ^ t* — 0 and for x £ 12. Therefore, we can continue 
the function q to the domain 12 x [ 0 , ti], where ti > t* and the properties of 
smoothness are preserved, and moreover, 

= * = 1,2, t<e. 

We shall construct the global temperature 9 in the form 

9 = e{x)f + g + (24) 

where q is the desired function, e(x) £ C^{[li,l 2 ]), and e = 1 for x £ 
[(,3i(0),(,ai(0)]. 


3 Construction of the weak asymptotic so¬ 
lution 


We fist consider the heat equation in system O- Using the formulas of the 
weak asymptotics method (see Section 5) and taking (|TT|) into account, we 
obtain 


du 

'm ^ 

+ 


^ , f^T{^20t - Viot) 

^(2 - B^o) +- ^ -2I00 

*^2*^0 D ^ , Pt {1^201 - ‘f lot) ^z 
-^(2 - ^ 60 ) +- ^2 - 


00 


(5(x - 

6{x 


^ 1 ) 

- (p2) + C>x)/(e), 


(25) 
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where the estimate is uniform in r, 

Boo= uJo{z)(jJo{-r] - z)dz, zu:o{z)uJo{-r] - z) dz, 

Jm Jm 

V = P(d, 

and the function p{t) is determined by the formula 

p{r) = (26) 

We assume that function q and its continuation are chosen so that the 
function q satisfies the boundary conditions of the original problem. Then the 
boundary conditions for the function q are zero. 

We substitute the function 6 determined by relation (24) and expression 
(25) for du/dt into Eq. (5). With accuracy Ox>'(e), we obtain 

- du 

L9 + — = LeT + Lq + Ai5{x - ipi) + A25{x - (/? 2 ), 

where 

q = i + l A = (-l)‘+‘^(2-Boo) + ^B5o. i = l,2. 

We let r tend to oo (i.e., for t <t*) and, in view of (10) and (18), obtain 

L(ef + g) = 2ip2t6{x - (^2) - ‘2(pitd{x - ipi). 

Since eT + q ^ 9 as t ^ 00 , we see that g ^ 0 as r ^ 00 . The total 
equation for the function q = q + q has the form 

L{q + e{x)I) = F{x, t, r) (27) 

- ^ ( 7 - Birm. - - if (. - ^.)i 

——']a - - ifi2) - H{x - .,31)], X G fi, 

where F{x, t, r) is a piecewise continuous function containing of terms that 
are uniformly bounded in e on x [0, ti] (they appear from Fi{x, t, r) in (21) 
and the derivatives of the function q). 

Equation (27) was derived under the assumption that (see (22)) 

2 

+T"] - A}d{x - ifi) = Ov'{e). (28) 

i=l 


(A^I- 
dx^ V 
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(29) 


By Qi, q* we denote the typical terms in q such that 

- ^p 2 ) - H{x - ipi)], 

Ip 

Lh = “ B)[H{x - ip 2 ) - H{x - (pi)], 

Lql = 2{x - x*)(7i+ - 72")V’“^^(r)[i^(x - P 2 ) - H{x - pi)], 
Lq 2 = -2(x - x*)( 7 f - 7^)V’“^-B(r)[i7(x - P 2 ) - H{x - (^ 1 )]. 


With accuracy up to functions smooth in 17, we can calculate them us¬ 
ing the fundamental solution of the heat equation (in what follows, we shall 
consider only the first equation), 


= 


. /•* (7i+(t0+72~(O) 

2 V^ Jo 

exp{-{x - ^f/{t -t')} 


r‘P'^ 

J 


Ip 


dt'di. 


(30) 


It is clear that the functions qi are uniformly bounded. It is also clear that 
qi 0 (7^(17), but qi G x [0,7i]) \ ({x = pi} U {x — <^ 2 })- Therefore, to 

justify Eq. (27), we must verify that no d-functions arise in calculating the 


derivatives 


d'^qi ^ ^ 

’ 'dx^ ^ (9t ’ dt ■ 

For this, we note that, for any test function C(x) up to functions smooth 

in 17, omitting the number factor, we have the relation 


{qi,C{x)) 



r*9?2 


V’ fitt 




t') dt] d^, 


(31) 


where C(j^) ~ solution of the heat equation 

ft — = C(C) at the point t — t'. It is clear that G C'°°(17) for all t. 

Calculating the derivative d‘^qi/dx‘^, we obtain 





V’ Vi(^,i 


t') dt' df,, 


where fi{f.,t) is a solution of the equation ft — f^^ = C'if,)- H is clear that 
the relation 

{§^xo))^Y.°ia<p.) + o{E) 

cannot hold for any coefficients Gj. The same is true for dqt/dt, dq*/dt. 
Similarly to (31), with accuracy up to smooth functions, we can write the 
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functions q*, i = 1,2. They have the same properties as qi with the only 
additional condition 

qi\x=x* = 0 . 

This easily follows from the explicit formula of the type of (31) and the fact 
that, for X = X* , the integral with respect to ^ is an integral of an odd 
function over a symmetric interval. Therefore, condition (28) is necessary for 
deriving Eq. (27). This implies that the function q + el is uniformly bonnded 
in X [0, ti] and belongs to for x ^ ipi, i = 1, 2. 

So if the function q + el satisfies Eq. (27) with zero initial conditions 
and the boundary conditions that follow from the fact that the function 9 = 
eT + q + q must satisfy the boundary conditions of the original problem, while 
the functions T and q are known (more precisely, they will be determined 
after the functions ipi are found), then we have the following assertion. 

Theorem 1. Suppose that the function 9 is determined by relation (28), and 
the function q is a solution of Eq. (27) with zero initial condition and zero 
boundary conditions. Suppose that relation (28) holds. 

Then the pair of functions 9 and u is a weak asymptotic solution of the 
heat equation in the phase field system, i.e., relation (5) holds. 

Remark 1. Expression (28) may seem to be strange, because we did not 
take into account the well-known fact that the Dirac functions are linearly 
independent. But this fact is taken into account in the framework of the weak 
asymptotics method (see Lemma 2, Section 5). Here we do not want to mix 
the substitution of the ansatz into the equation and the analysis of the results 
of this substitution (see Section 4). 


Now we consider the second equation of phase field system, i.e., the Allen- 
Cahn equation. By Definition [2l (see formula we must calculate the weak 
asymptotics of the expression 


„ def du 

dx 


eLu - _ ^9 


where the function 9 is determined by the relation (20). Calculating similarly 
as in Example 2 above (see Section 5 for details), we obtain 

jr =V^5{x-^i) + Vi5{x-^2) (32) 

+ Vi5'{x — if i) -\- V2d'{x — ^p 2 ) -\- Oxi'is), 


where Vt, i,j = 1,2, are linear combinations of several convolutions. Their 
expressions will be given below. Here we note the following. It is rather clear 
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that the expression 


£ 

2 


{uxfCx dx 


leads in the singular part to functions of the form 5'{x — (pi). 

Indeed, according to our construction, the expression eux is approximately 
the function (sign^z)^, where sign^ is the regularization of the function sign, 
i.e., a soliton-type function. It does not change its structure when squared, 
while its derivative divided by e becomes the derivative of d, i.e., 5'. 

Similarly, the term 



dx 


is transformed as 
f du 


dx 


{u — ■u'*)C dx = 


A 

dx 


F{u)( dx = — 


[ F{u)Cx 
Jr! 


dx. 


Similar arguments show that the singular expression itself in the asymptotics 
of the last integral is also of the form 6'. 

We denote 

^ {1 + ^oiz) + uJo{-z -rf)- uJo{z)uJo{-z - rj)} . (33) 

We have the following estimates: 


n{z,r]) = 1 + f{z,r])e‘^^, r] ^-oo, (34) 

n{z,ri) = uJo{z) + fi{z,ri)e~^'^, rj ^ oo, (35) 


where 


J \f (z, ri)\dz < const, 


j |/i(z, ry)| dz < const. 


These relations readily follow from the explicit form of the function ujq{z), 


. . e" - e " 

^^o(^) = ,, , 

+ e ^ 

In fact, relations (34) and (35) express the above-described properties of the 
ansatz u (11) in different terms. 

Using the technique of the weak asymptotics method (see Lemma 7 in 
Section 5), we can write the weak asymptotics of the expression Jj^i ^utUxdx 
in the form 


CutF 


2 

.dx = '^ (fit 
i=\ 




, 1 

^2 /32 



n'^)dz + Oj,,{e). 


(36) 


16 



Similarly, we obtain (see Lemmas 5 and 7, Section 5): 


/ 


duxC, dx = ]-[ipu+ {q + 


\x=ipi 


1 


- 2 \-V2t+ (g + g)L=, 


X=V?2 


^r)i^,v)C{^2)dz+ 0{£>^), /re (0,1/2). 

(37) 


Thus, adding expressions (36) and (37), we see that the coefficients of the 
(5-functions in formula (32) have the form 


= - 
^ i 

where 


^itBn + 


+ (el + q)\x=ipi)Cn, i —1,2, (38) 


B^ = 

[ n',iQ',-n'^)dz, 


Bh = ^ 

j[z{n'^ - 11(,) -iz + r/)^](f2/ - n'^) dz, 

(39) 

Cn = 

ji^'z - dz, q = q + q- 



Similarly, we obtain (see Lemma 6, Section 5): 


where 


vy = Vi 


(40) 


f 

/Ri 

(41) 

D = l 

[ F{n) dz 

/Rl 

(42) 


Thus, to obtain (3 = (5{rf) we have the equation (obviously, this is a neces¬ 
sary condition for the relation T = Ox)/(e) to hold, see (El): 


C 


(43) 


According to (41) and (42), the right-hand side of the relation (42) is positive. 
So we have proved the following assertion. 
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Theorem 2. Suppose that the assumptions of Theorem 1 and (43) are satisfied 
and 

2 

^ Vi5{x - <pi) + V^Six - (p 2 ) = Oofie). (44) 

i=l 

Then the pair of functions 6, u is a weak asymptotic solution of the phase 
field system o in the sense of Definition 2. 

Thus, under the assumption that the classical solution of the phase field 
system exists (see (i) and (ii) above), relations (27), (28), (43), and (44) are 
sufficient conditions for constructing a weak asymptotic solution of system (1). 
In what follows, we prove that the relations mentioned above are equations 
for determining the functions j3 = /3(t) and (/?ji(r), i = 1,2. We present an 
algorithm for solving these equations. 


4 An analysis of the relations obtained 


We begin with relation (28) and, for a while, forget everything said about the 
notion of linear independence. 

Then, for this relation to hold, it suffices to have the two relations 


Bbt + ii\ = 

In view of (18), we have + if = 
and adding relations (45), we obtain 


i = l,2. (45) 

Taking this into account 


x=(pi{t) 


b'o {pr- \ [ 

= B{t) I 

Here we used the relation 


w{z)uJo{-Z -v) { Pr - 


dx 


+ 


x=(pi{t) 


dx 


X=tf2{t) 


dz (46) 

= 2Bfib- 


bt = T2t - Tit = iPoPr + ^> 0(^21 - Tll)t, 


which follows from the definition of the function p, see (25). 

We set the last term to be zero, since we show below that {ip 2 i — Tii) = 
0(|r|“^). In view of Lemma 3 in Section 5, in this case we have {t 21 — 
Tii)tfio = C>{e). Moreover, in view of the Stefan conditions and the choice of 
the continuation of the functions and {poi(t), we have 


li + I 2 +7i +72 


dd 

dx 


+ 




de 

dx 


X=if>2{t) 
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We transform the left-hand side of relation (02); 


I = V-Oi 


= V'ot 

1 f 


1 


Pr 


UJo(z)uJo{-Z -r]) {pr- ^z ] dz 


W+o/ ^^oiz)uio{-z - r])^zdz 
ajo(z)uJoi-z - T]) (^Pr + -r]) + 


Next, we change the variables in the last integral 

Pr 


Ik 

k 


[ {-z - p)loq{z)loq{-z - p)dz =-^ [ zuJo{-z - p)u}oiz)dz 
Jki P Jri 


and note that 


Finally, we get 


f3r 


1 d 


I = V’Oi 


1 


PT 


kk 

2/3 dp 


= wk 


{zii>o{z)uJo{-z -p) + zuo(-z - p)u;o(z)} dz 

1 

UJq{z)uJo{ — Z — p) dz . 


Since the function ojq/z) is odd, we see that the expression in braces in the 
hrst integral in the right-hand side is 

d 

-z— (1 - U)o{z)u}q{z + p)). 


Hence, integrating by parts, we obtain 

/3r 


I = TpOt 


Pt 


2/32 Ai 


(1 - uJo{z)uJo{z + p)) dz 


1 dp 


(joo/z + p)uJo{z) dz 


Or, finally. 


where 


2/3 dr Jri 

I = V’Oi 

B{p) = f (1 — coq/z + p)loo{z)) dz = 277tanh?7. 


I d ~ 


(47) 
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So, in view of (46), we obtain the equation for p\ 


(48) 

Or 

d 

— + r/tanh?])) = 2B. 

After integration, we obtain 

??(1 + tanhr/) = 2/3 f B{T')dT'. (49) 

J — OO 

Since B{t) ^ 1 as t ^ oo, we have ^ 1 as r ^ oo. 

Since B{t) = for any N as t ^ —oo, the integral in the left-hand 

of (49) converges and hence the right-hand side tends to the limit 


lim rj (1 + tanhr]) 

T^ — OO 


lim 2(3 [ B{t') dr' = 0. 

J_oo 


Moreover, in view of the inequality i? > 0, we have rj > 0 for r G which 
implies that 

rj ^ 0 as T ^ —oo. 

Here we took into account the inequality /? > 0, which follows from (43). 

Substituting/? = ((71)“^( 77 ))^/^ into (49), we obtain the following equation 
for the function 77 : 

77 ( 1 -|-tanh? 7 ) = (CD”^)^/^ f B dr'. (50) 

Jo 

Since the functions contained in (50) are monotone and the limits exist as 
r ^ ± 00 , this equation is obviously solvable. 

Next, we have 


Pi'r) 



(3 = const, 


— 00 . 


Moreover, as readily follows from the exponential rate of convergence to the 
limit of the functions C/D, 1 — B, and 1 -|-tanh 77 as r ^ ± 00 , the derivatives 
satisfy the estimate 

= C>(|r|-^), |r|^oo. (51) 
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Now we calculate the limits of the expressions 


Ji = + %) - Ai, i = l,2, 

as r ^ ±oo. As r ^ oo, we have 


p^oo (p~r), 
Therefore, we have 




lim Ji = 0 


in view of the Stefan conditions (9), and Ji = 0{t ^), i = 1, 2. 
As r ^ — oo, we proved that ^ 0 and hence 


ipi-ip 2 = 0{e), T ^-oo. 


(52) 


Theorem 3. Relation (50) is a sufficient condition for relation (28) to hold. 

Proof. In view of the corollary of Lemma 2, Section 5, the estimate for the 
coefficients of the (5-functions in (28) as r ^ oo and (52) are sufficient condi¬ 
tions for relation (28) to follow from (46). In turn, relation (50) follows from 
(46). □ 

Now we analyze the relations = 0 and (44). We use the explicit form 
of the function F{z) and the formula 


CUo 



+ e-^■ 


(53) 


As rj ^ oo, we have 


n{z,r]) =(jJo{z)+ 0{e (54) 

and hence (5 ^ 1 as p ^ oo. 

From relation (33) we easily obtain 

l^l^oo, 

^^(2:,??) = C>(e"l^'), 1^1 ^00, 

Cn= u;o{z){'i--^oi-z-ri))dz >0. 

JRi 
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Similarly to (54), using (45), we can verify the relations 

= C>(e"^l^l), Tj^oo, 

Bn = l + 0{e-^^^^), v^oo, (55) 

Cq = 4 + 0(6“^’'), 7] ^ oo, 

Bn = 0, Tj^oo. 

We rewrite the expressions for the coefficients in more detail; 

Fi' = + Cn) + ^Bl, - ql^^Cn, 

^2 = —‘f2t{Bn + Cq) + B^ + 

Form formulas (57), the definition of the functions (pi, and the kinetic over¬ 
cooling conditions (9), we obtain 

lim F/ = 0, 

T^OO 

because q\x=ipi = 0, while g = 0 for t < t*, i.e., because of the fact that the 
model function T -|- el is constructed so that conditions (8) be satisfied. Thus, 
in view of the corollary of Lemma 2, Section 5, for the relation 

ViS{x - pi) + ViS{x - p) = Ov'{e) 

to hold, it is sufficient to have Vi + = 0 or, in more detail, 

(pit + P2t) {Bq + Co) — — q\x=ip2)Cn = 0. (56) 

Thus, we have proved the following assertion. 

Theorem 4. The conditions of Theorem 1 and relations (49), (55), (56) are 
sufficient for the functions 6, u to be a weak asymptotic solution of the phase 
field system. 

We consider relation (56). It follows from the above that Co = C^{rj) 
decreases sufficiently fast as r] ^ oo, in any case | |r/|Co(? 7 )| ^ const. This 
fact, the estimate 

q\x=<pi - q\x=ip 2 = C>{\pi - P 2 \^), //G (0,1/2), t<t*. 

proved in Lemma 4 in Section 5, relation (52) implies that the obvious estimate 
llv^i - P 2 \^Cn\ < econst|7?Co| = C(e^) 
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for t G [0, fi]. 

For t <t* , by Lemma 6, we have 

= o(iV’oV'ir), 

where tpi = (/? 2 i ~ V^ii- Hence, by (15), and a statement similar to Lemma 3 
in Section 5 (see also (17)), we have 

|2(v?iot + <^20t) - q\x=^i +q\x=^2 \ = t < t*. 

We introduce a function V{t) G (7“ such that G ^(M^), V{—oo) = 0, 
H(oo) = 1. Then, in view of considerations similar to those used in Lemma 3, 
we can show that the following estimate holds; 

(Ho + Cq)(V{T){ ipiQt + v^iot) + {q\x=(P2 ~ q\x=ipi)Cq = 0{e^), /x G (0,1/2), 
for t <t*. 

This implies that the left-hand side of (32) is estimated as Ox)'(e^) if 

{Bq + Cq)[—V (v?iot + </^20x) + ‘fit + ‘P2t\ — 2 Bq = 0. (57) 

From this relation we obtain 

+ ‘P2i)) = {y - l)(v^iot + V^20t) + ^Co ' 

From this equation we determine the extensions of the functions ^pn and 
P 2 i- We note that our argument results in an equation that does not contain 
the temperature, namely, the assumption that the classical solution exists 
till t = t* is sufficient for constructing a global smooth approximation of the 
solution. 

Let us calculate the function Bq = Bfi{rj) in more detail. We have 

Bn{v) = J ^'AK-K)dz. 

In this integral, we make the change of variable z ^ —z — r/. Then we obtain 

o;^-(L!'- h;), (0/- h;) ^-h/. 

Hence we have 

Bniri) = J Q'^fl'^dz. 
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Finally, we obtain 


Bn = 11 + dz = 


1 

2 


/ 


iK?dz. 


Thus, for finite ry, the denominator (Bq + Cq) in the last term in the right- 
hand side of (58) does not vanish. Moreover, using the explicit form of the 
functions ujq{z), B^, Bq, and Cq, we can verify that 

2^=0(|„|-''). W^oo, 

where ^ 1 is an arbitrary number. 

Therefore, we have 


+ <^12 


V^lOt + ^20t 


{V (r') - 1) dr + 


^ r /3; Bj, 

T io + ■ 


(59) 


We note that 

Pot (^1 + - V2ll)]^ = PoPr, (60) 

where p = and the functions r/ and (3 are determined by Eqs. (50) 

and (43). 

The system of Eqs. (61) and (62) allows one to find the functions ipn^ 
i = 1,2, and thus completely determines the functions contained in the ansatz 
of the weak asymptotic solution of system (1). 

It is easy to see that the solutions thus constructed satisfy conditions (15). 
Next, using the explicit formulas for B^, B^i, and Cq, we can easily verify 
that B^\fj=Q = 0 (because coo{z) is odd) and {BQ + Cfi)\rf=o / 0, because Coq{z) 
is even. This means that 


{p>it + ‘P2t)\-n=o = ^- (61) 

We note that the right-hand side of (27) is a piecewise smooth function 
for IV'I > const > 0 and smooth for x ^ ipi, i = 1,2. As tjt ^ Q {t ^ t*), 
the right-hand side of (27) becomes proportional to 5{x — x*), where Xq = 
x*\p=Q, and the proportionality coefficient is negative. Thus, if x and t vary 
in the respective neighborhoods of the points Xq and t*, then the d-function 
with a negative coefficient appears and disappears in the right-hand side of 
(27), which results in a negative soliton-like jump of the temperature in a 
neighborhood of the point t = t*, x = Xq. 
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Calculation of the temperature jump. To prove the statement in the Intro¬ 
duction concerning the temperature jump, we must calculate the quantity 

[el + q + . 

Let us verify that q + q is a continuous function. This function is the sum 
of solutions of the heat equation with singular right-hand sides, but these 
singularities arise after the substitution of the continuous function eT (which, 
of course, is not a solution) into the left-hand side. It remains to prove that 
the solution differs from the function eT by a continuous function. For this, it 
suffices to verify that the right-hand sides, which arise after the substitution 
of eT into the equation, do not generate any singularities in the solution of 
the heat equation in addition to those contained in the function eT. 

Here we, in contrast to the preceding statements, use the fact that the 
singularities in the right-hand sides arise as the result of the substitution. 

By qi, i = 1, 2, we denote the terms in q+q corresponding to the right-hand 
sides: 

-cpi)-H{x- ip2)], 

^ ~ “ *^^2) - H{x - V 9 i)], 


The other terms in q are solutions of the heat equation with piecewise con¬ 
tinuous right-hand side and hence are continuous. 

The functions qi have a similar property if they are calculated up to 0{e). 
Indeed, we denote 

^ - X){X - ip2) 

n = 7 - - - 

if 

and represent, for example, the function qi in the form 


1 r* B{T{t',e)) ^ - 

2 \/^ Jo Vt — t' 9 ^^ 



Omitting the number factor and integrating by parts in the integral over 
we obtain 



d^dt' 


Bile 


dt'. 
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The last term is a solution of the heat equation with the piecewise contin¬ 
uous right-hand side 

and hence it is continuous uniformly w.r.t e > 0. The other terms are also 
continous functions uniformly w.r.t e > 0. 

It is easy to see that in the formula for qi there is no term containing the 
derivative because this term would be of order 0(e). 

Indeed 


OB 

'dt 


Il[H(x - ipi) - H(x - <f 2 )] 




since |n| < Cxj: for x G [ipi,i.p 2 ] and e = p. The derivative decreases 
faster than any power of |/9!^| < const. This implies that 



Hence we have 

[q + Q\\x=x*,t=t* = 0. 


Let us calculate the function el. We have 

gI\x=x* = 2^^^* ~ ^2t) = ~^Pr- 

It follows from (48) that 


Pt T ^ oo, 

Pr ^ 0, T ^ —OO 

and hence ^ 

[el]\x=x*, = lim ((Plot + (P20t)- 

t=t* 2 

5 Technique of the weak asymptotics 
method 


First, we recall the definition of regularization of the generalized function. 
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Definition 3. A family of functions f{x, e) smooth for e > 0 and satisfying 
the condition 

w — lim/(x, e) = f(x) 

£—^0 

is called the regularization of the generalized function/(x). 

We note that, by definition, the last relation can be rewritten as 


liin(/(x,e),C(a:)) = (/,C) 

for any test function (from now on, ( , ) denotes the action of a generalized 
function on a test function). 

Lemma 1. Let Tt = {x — (p{t) = 0}, a: £ M, where (p{t) is a smooth function, 
let uj{z) £ S (§ is the Schwartz space), and let (3 = I3{t) > 0. Then the 
following relation holds for any test function C{x): 

e ^ 

where = ff°^Lo(z)dz. 

Proof. The expression in the right-hand side can be written as 

- f uj(/3- — C{x) dx = f u;{z) dz + 0{s). 

Here we perform the change of variables z = f3{x — (p))e and apply the Taylor 
formula to the integrand at the point x = ip. By definition, the last integral 
is the action of the generalized function (5~^Ai^6{x — p) on the test function 

C- □ 

Proof of the formula in Example 1 in the Introduction. Let lv{z) £ 
C°°, io' £ S(]R^), lim^^+oo u;(z) = 0, and lim^^-oo = 1. We verify that 


LO 


X — Xq 


— H{x — xo) = Oxi’is), xq = const. 


(62) 


where H is the Heaviside function. By Definition^ we consider the expression 


/ 


UJ 


X — Xq 


-H(x- xo) 


((x) dx 


(63) 


^C{xo) [ [uj{z) - H{z)] dz + 0{e‘^). 
JRI 
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Here we performed the change of variables z = {x — xo)/s and applied the 
Taylor formula to the functions ^it the point x = xq. In view of our 
assumptions and the properties of the function H, the integral in the right- 
hand side of the last relation converges and hence the right-hand side of (P|) 
is of order 0{e). We thus obtain estimate 

Proof of the formula in Example 2 in the Introduction. Let uji{z) G 
(cvi)' G §(M^), lim^^+oo Wi( 2 ) = 0, and lim 2 ^_oo= 1, i = 1,2. We 
consider the integral 


Integrating by parts, we obtain 

= V.. (/I "" 

" L (^) (^) 


We perform the change of variables z = (x — x^je with i = 1 in the first 
integral and with i = 2 in the second integral and apply the Taylor formula to 
the function F{x) = f^^Ciy)dy at the points x = Xi, i = 1,2, respectively. 
Then we calculate the first and the second integral, we have 


f Ax 


Ax 


J = Bi H{x - xi) -I- B 2 H{x - X 2 ) + Ox>'(e), Ax = xi - X2, 

Bi{p)= / (:oi{z)u]i{-z - p)dz, B 2 {p)= / uj 2 {z)uji{z - p)dz. 

./Ri JRi 


To calculate the linear combination of generalized functions up to Ox>'{£^), 
we must improve the classical definition of linear independence. This improve¬ 
ment plays the key role in the study of problems with interaction of nonlinear 
waves. 

Indeed, let (/>i / <1)2 be independent of x. We consider the expression 


gi5{x - (j)i) + g25{x - (t) 2 ) = a > 0, (64) 

where the functions gi are independent of e. Clearly, for the last relation to 
be satisfied, it suffices to have 


gi = 0{e^), i = l,2, 
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or, with the properties of the functions gi taken into account, 

gi = 0, i = l,2. 

But, if we assume that the coefficients g^ depend on the parameter e, then 
the above estimates do not work. Namely, let us consider the following specific 
case of this dependence: 

gi = Ai +Si{/\(t)/£), i = l,2, (65) 

where Ai are independent of e, the functions Si{a) decrease sufficiently fast 
as |(t| ^ oo, and A(^ = (/>2 — (/'i- 

Lemma 2 (Linear independence of generalized functions). Suppose that the 
estimate 

|(TS'j((T)| ^ const, i = 1,2, —oo < a < +oo 

holds. Then, for a = 1, expression ((Ml) implies the relations 

Ai = 0, A2 = 0, 51 + ^2 = 0. (66) 

Proof. Using the Taylor formula in dMI) and taking (|65|) into account, we 
obtain 

+ S2C{f^2) = 5lC(</-l) + S2Cif^l) + S2i^2 - 

where 0 < p < 1 Since the function aS2{(7) is uniformly bounded in u G 
we obtain 


*S’2(A0/e)((/>2 - (pi) = {-aS2{a)}\^^^^/^ • e = 0(e). 

So expression (IMl) can be rewritten as 

AiCih) + A 2 C(<P 2 ) + (Si + S 2 )C(<Pi) = 0(e). 

Thus, because the coefficients Ai are independent of e, we obtain the assertion 
of Lemma 12 □ 

Corollary 1. Suppose that 

|o'5'i(ci')| < const, i = 1,2, cr > 0 

and the functions (pi = <pi(t,£) are continuous and uniformly continuous and 
satisfy the condition that Acp > 0 for t > 0 and e > 0 and Acp = 0(e) for 
t < 0. Then, for a = 1, expression (66) implies (64). 
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Proof. Following the above argument, we obtain 

SiC(0i) + S2Cih) = (Si + S2)Ci^i) + S2ih - 0i)C'(0i + 

In view of our assumptions, 


S 2 {(f 2 - = 0{e) 


uniformly in t, which implies the desired assertion. □ 

Lemma 3. Suppose that f{t) £ C^, /(^o) = 0; f'{to) / 0. Suppose also 
that g{t, t) locally uniformly in t satisfies the condition 

\Tg{t,T)\ ^ const, \Tg^{t,T)\ ^ const, —oo < r < oo, 

and g{to,T) = 0. Then the inequality 


\ 9 {tJ{t)/e)\ ^ eQ, 


where = const, holds in any interval 0 ^ ^ t that does not contain zeros 

of the function f(t) except for Iq. 

Proof. The fraction is locally bounded in t. The fraction is also 

locally bounded. We have 




9 {tJ{t)/e) 
(t - to) 


f{t) _ t-tp 

e fit) ’ 


According to the assumptions of the lemma, the last factor in the right-hand 
side is bounded on the interval under study. The product of the hrst and 
second factors (without e) is bounded in view of the properties of the function 
g{t,T). □ 

Corollary 2. Suppose that the assumptions of Lemma 01 are satisfied for 
0 ^ T < oo (—oo < T ^ 0). Then the assertion of Lemma ^ holds on any 
half-interval which does not contain zeros of the function fit), and 

signt = sign/(t), t £ {to, t]. 

The proof of Corollary |2l is obvious. 


Lemma 4. The following inequality holds: 


9\x=ipi 9\x=ip2 — T2\^)- 
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Proof. We consider only one of the the functions q from (29), namely, the 
functions determined by relation (30), and apply the method developed in 
El We consider teh difference of expressions (30) omitting the number 
factors: 


\x=ipi q\x=ip2 


fl 

Jo Jip 




A{t - t') 


■ dr 


V31 - T 

I (<^2-<y?i)(</5i+ <y?2-20 

i m exp(-^^^^-^)[exp--1] 


4(t - P) 


4(t - t') 


tPl+y2 

2 


— t' 


^ ! (<^ l - 0 %r (<^1 - V ^ 2)(<^1 + <^2 -0 11 

t exp(-^y^— -1] 



4(t-d) 


4(t - 1') 


to V’Vt - 1' 

{t-t'Y {‘.p2-<^iY 


{p2-piY {t-t'Y 

where we choose /U G (0,1/2). 


■ dt', 


Next, using the inequality 


- 1 




< const 


for a > 0, X G [0, oo), and taking into account the fact that the integral 

{t - dt', /iG (0,1/2), 

converges, we obtain the statement of the lemma for ip 2 > If ^2 < 
(this is true for t > t*), then we must change the exponents outside the square 
brackets in the integrands. □ 



Lemma 5. Suppose that uj{z) G decreases faster than any power of jz 
as \z\ —> 00 . Then 


e ^u{{x - ipi)/e)qj = qj 6 {x - ipi) J u;{z) dz + Ov'{£^), 

where qj is one of the functions defined by the relations (29), i = 1,2, // G 

( 0 , 1 / 2 ). 
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Proof. We consider one of the functions q, namely, the function defined by 
the relations 


Q = 


1 dt' 


4{t -1') 


)df. 


The desired relation can be written as 


--1 


C(x)(j 


x-^Pi 

e 


qdx = C{Pi)q\x=^pi / ijj{z)dz + 0{£>") 


or, omitting the number factors, in the form 

{X - ifi)' 




C{x)uj 

dt' 


e 

r‘P2 


qdx 


(67) 


= C{pi) [ 

Jo 


0 - 1 ' 
* dt' 


exp 




rip2 


0 V'Vt - 1 ' 


exp 


VI 


{x-f,f 

4(t-f0 

4(t - t') 


df^ dx 

di + o{£''). 


In the left-hand side, we change the variables (x — (pi)/£ = -z and obtain 


d = J Ci‘Pi+£z)uj{z) j 


dt' 


rV2 


0 ikVt - t' 


exp 


VI 


4(t - t') 


dz df^. 


It is clear that to prove relation (67), it suffices to prove that, with accuracy 
up to small values, the term ez can be omitted in the exponent. 

We consider the expression 

'■* dt' 


J = y CiPi +£z)uj{z) j 
rv 

""J 


1py/t — t' 
4(t - t') 


exp 


A{t -1') 


dz df,. 


Since the function lv decreases fast, we can assume that \z\ < e~^, 6 G (0,1). 
Next it suffices to consider the integral over t' from 0 to f — because if 
t' G [t — then the obtained integral can be estimated as 0(6^^“^^/^). 

In the remaining integrals, we can use the same method as in the proof of 
Lemma 4. Namely, we transform the difference of the exponential functions; 

((/9i-C-Fe2;)2\ / _^)2' 


exp 


4(t - t') 


exp 


4(t - t') 


exp 
7 e (0,1/2), 


( 2 ( y’i - 0 - £z)£z 
4{t -1') 


- 1 


t — t' 


£Z 

t - t' 


-7 
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and apply the estimate 


^ exp ^ 


{2{ipi - 0 - ez)£z 
4(t - t') 


- 1 


t-t' 


ez 


/\ 7 


< const, 


which, obviously, is valid in the required range of variables. Hence, as in 
Lemma 4, we obtain the estimate 


J = 0{En, /xe (0,1/2), 

which proves relation (67). □ 

Similarly to Lemmas 4 and 5, we prove the following statement. 

Lemma 6. The following relation holds for any of the funtions qi defined by 
relations (29): 

qi{(pi,t) -qi{qio,t) = C’dV’oV'iD, h ^ (0,1/2). 

Lemma 7. Relations (38), (39) hold. 

Proof. In Section 3, we introduced the function Q.{z,rj). Of course, we can 
change the signs in the arguments using the fact that the function lvq{z) is 
odd, but the form presented above is convenient for calculations. 

Calculating the derivatives ut and Ux, we obtain 

lit = ^ { - x)]t cho + [/^(^ “ ^0 

- [I3{ti - x)]td}o ujo (68) 

- [I3{x - T2)]t‘^o } ’ 


Ux 


1 

2e 




+ djQ 




X-(p2 


(69) 


+wo 



/3 


X - (p 2 

e 


- dJo 


P 


X-(p2 


iOo 



Now we group the terms in the subintegral expression so that one group 
contain terms with the factor ujo{P{ipi — x)/e), and the other group contain 
terms with the factor Coq{P{x — ip 2 )/£)■ Thus, we write the integral as the sum 
of two integrals. Then we perform the change of variables 

x^Ti-j (70) 
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( 71 ) 


in the first integral and the change of variables 


ez 


in the second integral. 

From ( 68 ) and (69), we obtain 


utUxQ{x)dx 


(72) 


, /5r ,/ 


Pt 


^2t + ^V’O 


/32 


C((/?i)[(l + z){9.'^ -n'^) + {l-z- r])9.'^][Vl'^ - dz 
C(<7’2)[(1 + z)Vl'^^ + {l-z- dz 


+ 0{s). 


Here we used the formula for the derivative [3t = e~^f5r‘4’ot the fact 
that the functions in the integrand in (72) are invariant under the change 
z ^ —z — 77 , since each bracket is multiplied by —1 in this change. 

Now we calculate the other expressions contained in ( 6 ). We have 


J C'ix){uxfdx = ^ J 


-(Xo{ (3 


,7^1 - a; 


e 


cho /3 


,x - IP2 


+ Wo 

= f(C'(7?i)+ C'(7’2)) / {c:;o(2:)(l-wo(-z-7?)) 


,X -(P2 
e 


- Wo /3 


,X-ip2 


Wo (3 


e 

- X 

e 


- uJo{z)(Xo{-z wo(2;))(l - wo(- 2 : - rj))} dz + 0{e). 


(73) 

2 

dx 


It is easy to see that the integrand in the right-hand side of (73) can be written 
as 


J {(wo(2:)(l - Wo(-z - v)f 

- uJo{z)d}i-z wo(2:))(l - wo(- 2 ; - r/))} dz = ^ j{Q!^f dz. 

Thus, we finally obtain 

I J C'{u^?dx = j{C'{(pi) + CX^2)) J {n',fdz + 0{e). 

Then the left-hand side of (73) is positive, and hence the integrand expression 
in the right-hand side is also positive (until the measure of the points at which 
> const is positive). 
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The last term in (6) can be considered similarlly: 


1 

e 


F{u)Cx dx 


+C'iv>2)) 


F{Q)dz + 0{e). 


To prove this relation, we must successively perform the change of variables 
(71) in the integral in the left-hand side and then consider the half-sum of the 
integrals obtained, which, obviously, is equal to the original integral in the 
left-hand side. Moreover, we must take into account the estimates (34), (35), 
and 


n{z,r]) 


1-|-0(exp(—2^)), z ^ oo, 

ujo{—z — r]) + 0{ex.p{2z)), z^—oo, 


which together with the explicit form of the function F{u) 
imply that the integrals 


zF{Q) dz 



converge. 


□ 
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